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Abstract 



We study the role of direct (i.e. small-scale) instantons in QCD correlation 

functions for the nucleon. They generate sizeable, nonperturbative corrections 

to the conventional operator product expansion, which improve the quality of 

both QCD nucleon sum rules and cure the long-standing stability problem, 

in particular, of the chirally odd sum-rule. 
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QCD sum rules, introduced by Shifman, Vainshtein and Zakharov (SVZ) 0, provide a 
systematic, nonperturbative framework for the calculation of hadron properties. They have 
been intensely studied and applied over the last decade and produced the most exhaustive 
model-independent analysis of hadron properties to date. 

The sum-rule approach is based on the comparison of two "dual" descriptions for cor- 
relation functions of hadronic currents, in terms of quarks and gluons on the left-hand 
side (LHS) and in terms of hadrons on the right-hand side (RHS). The RHS uses a simple 
parametrization of the spectral function in terms of the hadron parameters, such as mass, 
overlap with the source current and continuum threshold. The QCD calculation on the LHS 
employs a non-perturbative operator-product expansion (OPE). Long-distance bulk proper- 
ties of the physical vacuum are efficiently parametrized in terms of the vacuum expectation 
values of composite quark-gluon operators ("condensates"), which are independent of the 
hadron considered. The short- distance physics is contained in the perturbatively calculated 
Wilson coefficients. The inverse renormalization scale of the operators, serves as the 
dividing line separating long and short distances. 

Both sides of the sum rules are then Borel-transformed, and the hadron parameters are 
determined by fitting the two sides in the fiducial region, i.e. in the range of Borel mass 
values in which both descriptions of the correlator are expected to be adequate. The quality 
of this fit is the only intrinsic criterion for the accuracy and reliability of the sum rules. If it 
is met sufficiently well, the resulting hadron parameters will be approximately independent 
of the Borel mass in the fiducial region. The nucleon sum rules, however, do not show such 
a stability plateau, despite many improvements over the last decade 0. It seems that some 
relevant physics in the fiducial region (around 1 GeV) is missing in the OPE. In this letter 
we suggest that small-size instantons 0, termed "direct" by SVZ, provide the dominant 
part of this physics. 

Instantons |1J are classical solutions of the euclidean Yang-Mills equation. Due to the 
infrared complexities of QCD, their quantum properties and vacuum distribution cannot 
yet be derived from first principles. A consistent picture of their importance and bulk 
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features has been established, however, by extensive phenomenological f5|||, analytical |7j 
and numerical studies (in the instanton liquid model and on the lattice, e.g. in ref. ||). 
They indicate, in particular, that the average instanton size p c in the vacuum is considerably 
smaller than the average separation R between instantons ||: 

p c ~ - fm, i?~lfm. (1) 

Most of the contributions of these rather small instantons to the correlation function are 
ignored in the conventional OPE: As the scale of these fields is smaller than the inverse 
renormalization scale (taken around /i ~ 0.5GeV), they would contribute to the Wilson 
coefficients, but do not show up in the perturbative evaluation. 

The aim of our paper is to calculate the leading instanton contributions to the nucleon 
correlator, using an instanton size-distribution || in accordance with the bulk features (|l]), 
and to study their effects in the nucleon sum rules. We do not use the detailed assumptions 
of the above-mentioned model calculations. 

The instanton contributions are mediated mainly by the quark zero-modes || in the in- 
stanton background field. Due to their particular chiral and color properties, the magnitude 
of instanton effects is channel dependent. Even if they can e.g. be safely neglected in the 
vector and axial-vector channel, they play a dominant role in the pseudoscalar sum rules 
|T0| (and more generally in the spin-0 channel), where the conventional OPE was known to 
fail 0. 

Our expectation of sizeable instanton effects in the nucleon sum rules is mainly based on 
parallels with the pseudoscalar sector. The interpolating fields (see below) in the nucleon 
correlator contain spin-0 diquarks, which receive zero-mode contributions of the same order 
as the pseudoscalars. Furthermore, the magnitude of the nucleon correlator at distances 
around 1 - 2 fm is much larger than the perturbative contribution JTTJ] , which is reminiscent 
of the strongly attractive correlations due to instantons in the pseudoscalar channel. 

Recently, Dorokhov and Kochelev |T2[ made a first attempt to calculate instanton con- 



tributions to the nucelon sum rule. However, as we will outline in the course of this paper, 



we do not agree with their results and many of their conclusions. 
The QCD nucleon sum-rules are based on the correlation function 

i J d 4 xe iqx <0\T V (x)rj(0)\0> = ill q {q 2 ) + U m (q 2 ), (2) 

evaluated at intermediate space-like momentum transfer, Q 2 = —q 2 ~ lGeV. We will use 
the most general interpolating field of the nucleon (with minimal mass dimension): 

7](x) = e abc [ [ua(x)Cd b (x)}"f 5 u c (x) + t [ul(x)C^f 5 d b (x)]u c (x) } . (3) 

Here, u, d denote the up and down quarks, C is the charge-conjugation Dirac matrix, and 
the real coefficient t specifies the linear combination of the currents containing scalar and 
pseudoscalar diquarks. Ioffe's current || corresponds to t = — 1. 

We intend to evaluate the dominant instanton contributions to the correlator eq. (Q) at 
distances relevant for the sum rules, g _1 ~ 0.2 fm jT|. As these distances are small compared 
to the average separation R ~ 1 fm between instantons, contributions from a single instanton 
should dominate multi-instanton effects. The small average size of the instantons further 
implies a sufficiently small gauge coupling g(Ap c ), which allows us to calculate the correlator 
in semiclassical approximation. 

To this end, we evaluate eq. (0) with the help of the quark propagator in an instanton 
background field ||, which is dominated by the zero-mode contribution 

± ) = ^iM^!M + 0(( rl) (4) 

m*(p) 

The superscript ± refers to an instanton/anti- instanton of fixed size p and position xq. The 
zero-mode solutions ip Q have been derived by 't Hooft in ref. |J. Due to interactions with 
QCD vacuum fields, the quarks acquire an effective mass, m*(p) = m q — \^ 2 p 2 < 0| qq |0 >, 
which has been given by SVZ (In the following, we will neglect the small current quark 
masses m q .) 

Due to the symmetry properties of the zero-mode solution, maximally two of the three 
quarks created by the interpolating field propagate in zero-mode states. We approximate the 
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propagation of the third quark in the continuum states [0] by the standard quark propagator 
used in the OPE and obtain 

< 0|Tr7(x)T7(0) |0 > = \ + ^ < 0| qq |0 >) 



X 



J dpn(p) 



P 1 



m* 2 (p) 



d 4 x 



((x - X ) 2 + p 2 ) 3 (Xo + p 2 ) 3 ' 



(5) 



where we define 



d = 6(t 2 - 1), c 2 = I \l3(t 2 + 1) + lOt 



(6) 



Due to translational and gauge invariance, the integration over all collective coordinates 
of the instanton except the scale p is straightforward. The broken scale invariance, however, 
gives rise to a p-dependent weighting factor, the instanton size-density n(p). As motivated 
above, we will use the simple form proposed in ref. 0], 



n(p) = n c 5(p- p c ), 



(7) 



which incorporates the already described features and scales of the instanton distribution 
and the quite sharply peaked (half-width ~ O.lfm), almost gaussian result of Monte-Carlo 



simulations [14 



The Fourier transform of (|5|) contains the instanton contribution to the invariant am- 
plitudes H g (q 2 ) and II m (g 2 ) in (0). Their subsequent Borel transform, required in the sum 
rules, leads to quite complicated integrals ]TB| which we evaluate numerically, but do not 
write down explicitly here. Instead we present a saddle-point approximation, which is more 
transparent and allows for a direct comparison with the work of 
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where Mb is the Borel mass, erfc(z) is the complementary error function |Tg and we define 
z = MbPc- Note that the instanton contribution to Tl q vanishes for the Ioffe current (i.e. 
ci = 0). 



We eliminate the n c dependence of (|8|) and (0) with the self-consistency relation JT7| 



f nip) n r 
< qq >= -2 / dp -%L = 10 
J m*(p) m*(p c ) 

While the assumption of a completely zero-mode-induced quark condensate is rather strong, 



our final results will not depend on it, because the numerical value obtained from ( flOf) is 
almost exactly equal to the phenomenologically estimated value. 

Before combining eqs. (§) and (§) with the standard OPE terms of the nucleon sum rules 
0, we have to address the issue of double counting instanton physics. The conventional OPE 
already accounts for all nonpertubative vacuum fields (including instantons) with scales below 
the renormalization scale p in the condensates, whereas the whole M-dependence, logarithms 
and powers of M~ 2 , comes from the perturbatively calculated Wilson coefficients. 

The M-dependence of the instanton contributions, however, is exponential in all terms 
of Il m and in all but two terms of Il q . Hence these terms constitute new, nonpertubative 
contributions to the Wilson coefficients, originating from small- instanton physics. Il q , how- 
ever, contains also two power terms ||18|| , which represent instanton contributions to the 
condensates of dimension 6 (four-quark condensates) and 8 [ljj. In this letter we will re- 
strict the conventional OPE to operators of dimension < 6, so that double counting has to 
be prevented only for the four-quark condensates. 

The standard sum rule practice is to approximate the four-quark condensates by fac- 
torizing them into two-quark condensates, and it is known that this approximation doesn't 
always work well for the instanton contribution. Indeed, the four-quark condensate term in 
eq. d§|) is considerably larger than the corresponding factorized OPE condensate, which we 
will therefore omit completely. 

The LHS of the sum rules can now be obtained by adding the invariant amplitudes, (g) 
and (H), to their standard OPE counterpart with the corrected four-quark condensate. 
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On the RHS, these amplitudes are written in a (Borel-transformed) Lehmann representation 
with the usual pole-continuum approximation for the spectral functions ||. Equating both 
sides and transferring the continuum contribution to the left, we obtain the two sum rules 



Q 
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2c 3 M b A 2 (M, W) + 2c 3 7r 2 M 2 < -G^GT > A (M, W) + -c 5 (4tt) 4 < qq > 2 +—c lP , 
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and 



- (An) 1 <qq> 



c 4 M 4 Ai(M, WO - |-m^M 2 A (M, + csv^Pe" 4 ^ (z 2 - | j c 



mA^e-™ 2 ^, (12; 



where the 

A n (M, W) = \- e - w2 / M2 
contain the continuum contributions and 



1 vi^r' 



(13) 
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Let us now compare our results with those of ref. [p^ , where semi-classical instanton con- 
tributions to the nucleon sum rules are also evaluated. We find four additional terms in Ilg 
and one in Il m , all with exponential Borel-mass dependence. These terms, which apparently 
have been missed in ||12|| , are in no way subleading and their omission has drastic conse- 
quences: 1.) The spectral function of Il q becomes negative in the whole fiducial domain, 
thus violating its very general positivity bound |[20|| . 2.) For the same reason, the (usually 
more reliable) Ilq sum rule becomes unstable and has to be abandoned. Half of the physical 
information and the valuable consistency check between the two sum rules is lost. 3.) A 
numerically important part of tl m is missing and the relation between nucleon mass and 
quark condensate, usually manifest in the "Ioffe formula" 0, is obscured. 



Furthermore, the authors of |12| used the (only qualitatively reliable) saddle-point ap- 
proximation to evaluate the sum rules, which induces substantial errors (~ 30 %) in the 
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nucleoli parameters. (As already mentioned, our results are based on the exact, numer- 
ical evaluation of the instanton contribution.) Finally, double-counting of the four-quark 



condensate is not addressed in |12| and the 4- and 5-dimensional condensates are neglected. 

The free parameters to be determined in the sum rules are the continuum threshold W, 
the nucleon mass m and A at, the overlap of the wave packet created by r] with the nucleon 
(< 0| rj(0) \N a (k) >= \NU a (k), Xn = (47t) 2 Atv). We obtain their values by minimizing the 
difference, defined by the measure 5 of ref . @ , between the two sides of the sum rules (|TT1) 
and (|i~2l) in the fiducial Borel mass domain 0.8GeV < M < 1.2GeV. In this letter we choose 
t = —1.1 ||21|| , corresponding to the optimal current of Espriu et al. H and close to the 
Ioffe current. For any choice of the current, however, the instanton contribution will be 
substantial in at least one sum rule. 

Figure la shows the OPE contribution, the instanton contribution and their sum, the 
complete LHS, in comparison with the pole term A^e"™ 2 ^ 2 for the $ sum rule (|ll|). Figure 
lb shows the same curves, divided by the nucleon mass m, for the sum rule flT2"|). As 
expected, the instanton contributions to the ^ sum-rule are rather small for our current, but 
still improve the fit quality (S = 1.2%). The nucleon parameters become m = 0.88GeV, 
W = 1.43GeV and ~\% = 2.90GeV 6 . 

In the traditionally less stable m sum-rule, however, the instanton contributions become 
comparable to the OPE contributions and improve the fit quality dramatically. Their M 
behaviour balances the OPE terms over the full Borel-mass interval to bring their sum very 
close (5 = 0.03%) to the pole contribution. The nucleon parameters become m = 0.90GeV, 
W = 1.60GeV and A^ = 5.68GeV 6 . 

The different magnitude of the instanton contributions in the two sum rules implies that 
a simultaneous fit of both sum rules will not yield their individually optimal parameter 
values. Indeed, the nucleon mass increases (m = 1.06GeV) and the fit quality does not 
improve as significantly as in the m sum-rule. This might indicate some missing physics in 
the <ji sum-rule. Radiative 0(a s ) corrections to the conventional OPE, e.g., are known to 
contribute at the same order as the leading perturbative term to the ^ sum-rule and their 
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inclusion lowers the nucleon mass |22] 



As mentioned above, the quality of the sum rules is reflected in the approximate Borel- 
mass independence of m(M) and A^(M), obtained by solving the optimized sum rules. The 
instanton-improved sum rules leave both quantities quite insensitive to the Borel mass (Fig. 
3). The usually less stable m sum- rule, which receives the main instanton corrections, shows 
a particularly pronounced improvement and generates perfect stability plateaus over the 
whole fiducial region. 

The extension of our investigation to other baryons and to finite density and temperature 
is in progress. This work was supported in part by the US Department of Energy under 
Grant No. DE-FG02-93ER-40762. 
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FIGURES 

FIG. 1. a) The OPE (dashed line) and instanton (dot-dashed line) contributions to the $ sum 
rule. Their sum (dotted line) is compared to the pole contributions (solid line), b) The same for 
the m sum rule. The sum of OPE and instanton contributions is practically indistinguishable from 
the pole term. 

FIG. 2. The nucleon mass as a function of the Borel mass from (a) the $ sum rule, (b) the 
m sum-rule and (c) the combined fit. The coupling from (d) the $ sum-rule and (e) the m 
sum-rule. 
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